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1. Introduction 

The understanding of the structure of nucleon has been one of the central issues in hadron 
physics. For instance, the parton distribution functions (PDFs) have been studied extensively, and 
the observation of scaling violation in PDFs provides the cradle for the establishment of the fun¬ 
damental theory, QCD. Yet, there exist many unresolved questions for the structure of the nucleon. 
Ever since the EMC experiment showed that the proton spin carried by quarks is small[l], large 
effort has been made in both the experiment and the theoretical frontiers to identify all the contri¬ 
butions to the nucleon spin. Calculating the momentum fraction {x) is an integral part of the study 
of this subject. 

However, it has been found that calculating the disconnected insertion contribution to this 
quantity, which is necessary for the strange momentum fraction, is extremely difficult. Most lattice 
calculations are done for the connected insertion contribution. So far, only one quenched calcu¬ 
lation exists that includes all contributions to the nucleon spin[ 2 ]. In this paper, we calculate the 
strange momentum fraction with the overlap fermion on 2 -|- 1 flavor domain wall fermion configu¬ 
rations with the help of an array of lattice techniques. We also calculate the momentum fraction in 
the u/d channel for the disconnected insertion (DI), and take the ratio between strange and u/d chan¬ 
nel in DI. In the end, we compare our result with previous lattice calculations, as well as current 
global analyses of parton distribution at small x. 

2. Formalism and simulation parameters 

The momentum fraction carried by quarks can be obtained by calculating the following quark 
operator of the QCD energy-momentum tensor in the forward matrix element of the nucleon at 
finite momentum. 

Tii = - ^ (^qY 4 Di q + qfi Da (2.1) 

in which 



With this operator, we calculate the disconnected three-point function 


C3{tl,t2,p) = {PN{t2,p)L4i{tl))—C 2 it2,p){L4i{ti)), ti < t2 (2-2) 

in which the angular bracket denotes ensemble average, Pn is the nucleon propagator, and L denotes 
the energy-momentum tensor contracted in a quark loop. The momentum-projected nucleon two- 
point function is 


C^(i 2 ,p) = (Pv(f 2 ,p)) = £e-'>"(0|r[z(T,t2)z(0,0)]|0). (2.3) 

X 

where the nucleon interpolation field. Note that we take the quark source to be at the 

origin x = 0 , t = 0 in the above expression. In the actual calculation, we take the source to be a 
smeared grid source with a Z3 noise. A diagram representing the three-point function is drawn in 
figure 1 
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Figure 1: The disconnected insertion 


The momentum fraction can then be obtained via the following ratio 


1 Tr[reC3(ti,t2,p)] , V 

p Tr[r,C^{t2,p)] 


t2 > > 0 


(2.4) 


in which Fg is the parity projection operator (1 + 74 )/ 2 . To improve signal/noise ratio, we use the 
sum method[3] by summing the insertion time ti over an appropriate range 


ly Tr[r,C3(ti,t2,p)] 

Pt^[ T'A^eC^{t2,P)] 


—> (v)t 2 + const. 


(2.5) 


We use the overlap fermion valence on Nf = 2 + 1 domain wall fermion configurations by the 
RBC/UKQCD collaboration[4], for which the lattice spacing is = l.VVGeV. The light quark 
mass is 0.005 in lattice unit, which corresponds to nijt = 330MeV. For this work, 200 configurations 
are used. In the valence sector, quark masses iriq = 0.00809, 0.0135, 0.016, 0.0203, and 0.0576 are 
chosen to be included. 

On the choice of momentum for the two-point function, in principle, one may use any finite 
momentum allowed on the lattice. We choose momenta p = InlLa and An/La, and take an appro¬ 
priate average to acquire a final resulf. 


3. Lattice techniques 

For fhe fwo-poinf correlafion function, several fechniques are applied fo improve fhe sig- 
naVnoise rafio. We use an 8 -grid source wifh gaussian smearing[5] fo enhance fhe signal of fhe 
ground sfafe and Z 3 stochastic noise is used fo fie fhe fhree-quarks in each smeared source fogefher 
fo have a beffer overlap wifh fhe nucleon. Low mode subsfifulion[ 6 , 7] is also utilized; low mode 
confribufion fo fhe fwo-poinf funcfion is calculafed exacfly. Furthermore, for each configurafion, 
we use 32 sources, each of which is shiffed in fhe fime direcfion. Since af large t where fhe nucleon 
effeclive mass plafeau sfaits fo appear, high mode confribufion is less imporfanf, fhe inversion of 
fhe quark mafrix wifh low-mode deflalion is calculafed wifh lower precision to save inversion lime. 
Figure 2 shows fhe difference in eiTor bars of a poinl source, grid source, and grid source wifh low 
mode subslilulion (all gaussian-smeared). We should mention lhal wilhoul low-mode subslilulion, 
fhe signal from fhe grid source is even worse lhan lhal from fhe smeared poinl source, despite 
having more slalislics (also see [ 6 ]). 
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Figure 2: Comparison of nucleon effective mass at nijc = 330 MeV plotted from different source 
schemes (as ordered in the legend: a single point, 8 -grid without LMS, 8 -grid with LMS, and 8 -grid 
with LMS from 32 sources per configuration as described in this paper). In this graph, all sources 
are gaussian smeared. aAIJGrid w/ LMS & 32 srcs avgaAi is used in this work. 
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Figure 3: Source scheme for the loop. The right figure shows fhe grid poinfs on differenf time 
slices in each dilufed inversion. 


For fhe loop, we use a grid source wifh Z 4 noises in bofh spacial and time direction, and fhe 
sources are arranged in a even-odd fashion. For each configuration, 32 sources are used. Among 
all fhese sources, fhe grid is shifted in spacial direcfions, so fhaf fhese sources fogefher cover all 
time slices, and 1/4 of fhe lattice sites in any given lime slice. A skefch of fhis source scheme is 
drawn in figure 3. Similar fo fhe Iwo-poinl functions, we utilize fhe low mode averaging (LMA) 
technique[ 8 ], and calculate fhe low mode conlribulion lo fhe loop exacfly wilhouf any slochaslic 
noise, and high mode conlribulion is calculated wilh lower precision. 
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(a) {x) slope in the u/d channel (b) {x) slope in the strange channel 

Figure 4: 3-point/2-point function ratio in the form of slope as a function of the nucleon sink-source 
time separation from the sum method. 

4. Results 

In order to improve the signal/noise ratio, we use the sum method to calculate the physical 
quantities we want, i.e. {x), from the 3-point function (eq. 2.5). Figure 4a shows the signals in 
the u/d channel, in which both the 2-point function and the loop are with quark mass 0.016 (nij^ = 
330 MeV); figure 4b shows the signal in the strange channel, in which the nucleon propagator is 
calculated with quark mass 0.016, and the loop with 0.0576 which is at the strange quark mass. All 
the techniques described above contribute to the resulting good signals. In the u/d channel, we are 
able to reach a clear signal almost 9 sigmas away from zero. Furthermore, contribution from high 
modes and low modes are also plotted. It is evident that the high modes contribution is significant 
for {x ), hence we devote a large amount of computational time in order to control its error. 

We show in figure 5 and {x)s as a function of for several quark masses in fhe nucleon 

propagafor and fhe u/d quark loop; while fhe sfrange quark mass in fhe loop is fixed. We fhen 
exfrapolafe fo fhe physical pion mass using fhe following functional forms respecfively 

{x)^}=ai+a2m\ + a^m\, {x)s = bi+b2m\. (4.1) 

in which fhe a’s and /j’s are tiffing paramefers. Thus we obfain bofh and {x)s af fhe physical 
poinf 

= 0.0285(57), (x), = 0.0195(26). 

They will be renormalized and mafched fo fhe MS scheme [9]. 

In figure 6, we plof fhe ratio {x)s /for each quark mass, and fake a chiral exfrapolalion 
in ferms of fhe linear m|. The exfrapolafed value is found fo be 

='>-78P)- 

\Muld 
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Figure 5: {x)ud and {x)s extrapolated to the physical value 
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Figure 6 : {x)s/ ratio extrapolated to the physical pion mass. 
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5. Conclusion and Discussion 


We have studied the momentum fraction carried by quarks in nucleon with overlap fermions 
on 2 + 1 flavors domain wall dynamical fermion configurations. We utilize an array of lattice 
techniques including low mode substitution and low mode averaging, as well as grid source with 
Z 3 and Z 4 noises to greatly improve the signal of both the nucleon propagator and the quark loop. 
We find clear signals for {x)^j^ with a 5 sigma signal and {x)s with a 7 sigma signal at the physical 
pion mass. It is worth noting that our ratio result {x)s/{x)^j^ = 0.78(3) is consistent with the global 
analysis shown in figure 7 [10], given that the “disconnected sea” contribution dominates small x 
region due to its dependence[ll, 12 ], and that the {s + s)/{u + d) ratio is relatively flat in 
small X region. We would also like to point out that our result is in agreement with previous lattice 
calculations by T. Doi et al. [13] and M. Deka et al. [2]. This ratio can be used to better constrain 
the global fitting of the parton distribution functions at small x. 
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(s+s)/(u+d) distribution at = 5 GeV^ 



Figure 7 : Ratio of 5 + 5 over u + d versus x at = 5 GeV^ from various recent PDFs. 
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